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We have calculated the contribution of intersubband 
transitions to the third order optical nonlinear susceptibility, 
\^ (uj , u , u) for nonresonant as well as resonant third har- 
monic generation and \ (w, —lu,uj) for nonlinear refraction 
and absorption. As examples, we consider InAs/AlSb and 
GaAs/GaAlAs quantum wells. The effects of finite barrier 
height, energy band nonparabolicity, and high carrier concen- 
trations are included. It is shown that quantum confinement, 
rather than the band nonparabolicity, is responsible for high 
values of nonresonant Very high values of XzSzz are 

obtained for third harmonic generation and two photon ab- 
sorption for incident wavelength near 10.6 /im. Intensity de- 
pendence of refractive index and of absorption co-efficient is 
also discussed for intensity well above the saturation intensity. 
Effective medium theory is used to incorporate the collective 
effects. 



I. INTRODUCTION 

Intersubband transitions (ISBT) in semiconductor 
quantum wells (QW) have attracted a lot of attention!!!,. J 
in nonlinear optics since the work of West and Eglashi. 
As emphasized in a recent review by Almogy and Yarivu, 
the main reason for this interest is the occurrence of nar- 
row transitions with large oscillator strengths. In bulk 
semiconductors mobile carrier optical nonlinearity, x^ 
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due to intrabandtransitions is proportional to 
low frequenciesQ'El, where photon energy Uu) <C E g , the 
semiconductor bandgap; E is the electron energy, and k is 
the electron wave vector. This suggests that nonlinearity 
in bulk semiconductors increases with increase in the en- 
ergy band nonparabolicity. Some attemptsLrQ have been 
made to extend the above relation to calculate x^ due 
to mobile carriers in a semiconductor QW with a view to 
enhance the nonparabolicity and hence the optical no 
linearity. However, as noted by Jha and Bloemberge 
the low frequency limit is valid when tko is much smaller 
than the relevant bandgaps, which in the present case 
is the intersubband gap, Ejsbt- Thus in a superlatticc 
(SL) or a multiple quantum well (MQW) there are two 
relevant low-frequency regimes: Ejsbt < < E g , and 
second, Tilj < Ejsbt- For the former regime, which is 
of more practical interest, the above simple relation is 
not applicable and it is important to understand the role 



played by quantum confinement ! 1 °H 1J I Resonant enhance- 
ment of nonlinearity at ISBT is the second major theme 
in the exploitation of quantum wells as nonlinear materi- 
als. High values of x^ have been reported by thirdr-hp 
monic generation (THG) involving ISBT resonancestaE 
The third important aspect of optical nonlinearities in- 
volving ISBT is the significant effects of electron-electron 
interaction. Two of the most important consequences of 
including the Coulomb interaction are the shift of the-iit. 
tersubband resonance due to depolarization effects — Jill 
and the predicted occurence of optical bistabilityMlil 
Most of the work cited above has emphasized one or the 
other aspect. In this paper we develope a theory of the 
optical nonlinearities in the intersubband transition re- 
gion including the essential complexities such as, finite 
barrier height which allows the electron wavefunctions 
to extend into the barrier region, nonparabolicity which 
affects the intersubband energies, the subband effective 
masses, and the collective effect which is important at 
high carrier concentrations. To do this we first calculate 
the single particle response including realistic subband 
dispersions and then include electron-electron interaction 
through an effective medium theory. 

In Section II we present the theoretical formulation. 
Various approximations necessary to deal with non- 
parabolic bands in quantum wells are stated as well the 
necessary details of the expressions used for evaluating 
nonlinear susceptibilities. The nonlinear susceptibility 
X^zzzz{— 3w, lo, uj, uj), here after denoted as x^{3uj) de- 
scribing the THG and Xzzzz(— u>,u>, —uj,uj), hereafter de- 
noted as x^ 3 \— w ) are treated separately. To deal with 
resonant effects at high intensities power broadening ef- 
fects in resonances are also included. We also describe 
here the effective medium theory to account for depolar- 
ization terms due to electron-electron interaction. 

In Section III the results of numerical calculations of 
THG involving ISBT are discussed for InAs/AlSb MQW 
systems. This material system has become important 
because of its high barrier height (1.35 eV) which allows 
large charge density (of the order of 10 19 cm -3 ) to ac- 
cumulate. Very high room temperature mobility (2 x 
10 4 cm 2 V~~ 1 s~ 1 ), and ultra-fast energy relaxation time 
(1 ps at room temperature) of the system are also useful 
for device applicationsEj. We show a very high value of 
|X^(3w)| (of the order of 10~ 4 e.s.u. at 10.6 [mi) can be 
achieved from this system. 
In Section IV NRA, x (3) ( 

—a;), are calculated for 
the same material system and also for GaAs/GaAlAs 
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MQW system. We find that large value of x^ 3 H~ uj) 
can be achieved at two photon absorption (TPA) wave- 
length which is near ~ 10.6/zm for InAs/AlSb system. 
Since the values of Re(x < - 3 - > (— w)) is larger than that of 
lm(x^ {— ^)) near TPA wavelength useful applications 
of these systems in optical switching devices may be ep. 
pected. Optical bistability has already been predictedMI 
for InAs / AlSb system after considering electron-electron 
interaction through the effective medium theory. This in- 
trinsic optical bistability of a medium is of considerable 
theoretical interest. 

In Section V the numerical results of intensity- 
dependent absorption including the power-broadening 
term are presented. The calculated values obtained 
by using the effective medium theory to include the 
depolarization., effects compare, very well with earlier 
experimentalE!! and theoreticaO results. 

The conclusions are presented in Section VI. 



and 



m* Ei (E)=m*{0){l + ai(E-E ci )}, (3a) 



mUE) = m|(0) {1 + 2a t (E ~ E cl )} , (3b) 



where m*(0) is the band-edge effective mass in the bulk 
material, and on is the nonparabolicity parameter of the 
i layer. 

We calculate the third order optical susceptibility, 

(3) 

Xzzzz in the z-direction. In this paper susceptibilities 
in the in-plane directions will not be considered. The ex- 
pression for xfjli ( _w cr , wi, u)2, CJ3), where u a = w\ + lu 2 + 
U3 can be obtained by using the standard density matrix 
perturbation formalisms, and has 48 terms which are 
reduced to 24 distinct terms for i=j = k = l = z. 

Resonant third order susceptibilities have been calcu- 
lated for the cases of THG and NRA by the sum-over- 
states method. For THG we get, 



II. THEORETICAL FORMULATION OF x (3) 

In this section we first formulate the single particle 
response x^ f° r ISBT, and then the electron-electron 
interaction is included through the effective medium the- 
ory. 



Single particle response 

We consider a MQW system consisting of rectangular 
wells and barriers of widths L\y and L B , respectively. 
The well region is doped n-type so that at least the lowest 
subband is always populated. The z- direction is chosen 
to be perpendicular to the plane of the wells with z = 
at the center of a well. The ecpation for the envelope 
function, F(z) is then written asE3, 



d 2 Fi(z) 



2m* E .{E) dz 2 



(E - E ci ) 



n 2 kj 

2m* Ei (E) 



Fi{z) = 0, 
(1) 



and is solved with the continuity conditions, 

F w (z)=F B {z), (2a) 



and 



dFw (z) 



dF B {z) 



Til 



vw 



(E) dz 



m 



V B 



(E) dz 



(2b) 



for \z\ — Lw/2- In Eqs. (|l|) and (g) E ci represents the 
band-edge energy, E is the energy eigenvalue including 
the component due to in-plane motion, i stands for ei- 
ther the well (W) or the barrier (B), hk± is the in-plane 



momentum of the electron; 



and 



•V, 



arc the energy 



and the velocity effective masses, respectively, and are 
given by 
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where f(E) is the Fermi function given by 
f(E) ' 



1 + exp{(E -E a - E f )/k B TY 



(5) 



where E is the energy corresponding to the in-plane mo- 
mentum of the electron, E a is the energy eigenvalue of 
the state a of the QW, and Ef is the Fermi energy. The 
factor S(u) in Eq. (|) is 

S(u>) = E EX Z abZbcZcdZda 
b^a c^a d^a 

1 



(uj ba - 3uj)(u} ca - 2u)(uj da - lo) 
1 

( W ba + - 2u>)(uj da - lo) 

1 

Ka + W )( W c*a + 2w)(^rfa - w) 
1 

Ka+^)«a + 2w)Ka + 3w) 

EE I**IW 



1 



2UJ \ (u>da ~ u){LU ba - 3w) 
1 



(G) 



where Zjk are dipole matrix elements between two enve- 
lope states in the j-th and fc-th subbands. In Eq. (||), 
hujjk = Ejk —iTjk, where Ejk is the intersubband energy 
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between the states j and k, and is the broadening 
parameter; and oj* k is the complex conjugate of ojjk- In 
general, Zjk and Ejk depend on the in-plane momentum 
of the electron for nonparabolic bands. 

The NRA are the other important processes which 
has not been discussed in detail in the litcraturco. 
The importance of a complete theoretical expression of 
Xzzzz ( — WpW, —oj, oj) was pointed out earlier by Bloember- 
gen et. alHj for a system where close one and two photon 
resonances were present. In this article we present a com- 
plete expression of Xzzzz{~ oj,oj, —oj,oj) due to ISBT for 
a heavily doped MQW system in which resonances axa 
similar to the case considered by Bloembergen et. al.E£l 
The third order susceptibility for NRA is given by 

Xzzzz{-U,UJ,-UJ,Uj) -- 



where 



T(U>) = X! X! X! Z abZbcZcdZda 
b^a c^a d^a 
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The terms following the double summations in the ex- 
pressions for S(u>) and T(lo) are the seperated out terms 
containing the ground state as an intermediate state 
(c = a). The apparently divergent terms (c = a) in.-thc 
expression for T(oj) are reduced in the usual manncnHj. 

It should be noted that the expression for x*- 3 -^— oj) is 
complete and is rigorously calculated from the original 24 
terms which are obtained from the perturbation theory 
in dissipative systems. These 24 terms then reduce to 12 
terms due to io\ — —0J2 — ^3 = oj. Since the imaginary 
part of x^ 3 H~ w ) gives the two photon absorption (TPA) 
co-efficients, it should be positive for the entire frequency 
range. Taking the proper sign of the damping factor 
is, therefore, important. In the calculation of THG, the 
terms are calculated with the help of the perturbation 
theory in the absence of damping, and then the damping 
factors are introduced at the final expression to obtain 
the LorentzianJproadening at oj, 2oj, and 3oj resonancescJ. 
It was shownED that in the THG dispersion, the term 
containing [{oj ba — 3uj)(oj ca — 2oj)(ojd a — ^0] is dominant 
among all the 48 terms and no cancellation of resonant 
terms is involved. Hence our approximation for is 
justified. 

The nonresonant optical nonlinearity can be calculated 
from Eq. by putting r = taking the limit, oj — > 0. 
The resultant expression is 



Y (3) 

A zzzz 



(0) 



4e 4 



h ttL 



w Jo 



E 



m* Vw (E)f(E) 



3 



EEE 

bj^a c^a d^a 



Z-ab^bc ZcdZda 
^ba^ca^da 



EE 

b^a d^a 



\Zab\ 2 \z d a\ 2 
UJda^ba 



dE 



(9) 



In the infinite barrier height limit, the above expres- 
sion may be cross-checked with the formula derived by 
Rustagi and DucuingEa for the particle-in-a-box model, 
which when suitably modified for an infinitely deep semi- 
conductor QW takes the following form: 



p 4 N m* 3 T 9 JL 
(3) (q\ _ e ^sm w -L w \ -» 

Xzzzz \ ) 
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(10) 



where N s is the sheet carrier concentration, m^y is the 
bulk bandedge effective mass and N is the total number 
of occupied subbands. 

To calculate the intensity-dependent absorption we 
have included power broadening following the general 
formalisms to obtain the expression for the total (linear 
and nonlinear combined) optical susceptibility, x(— u>) as: 



the effective medium theory to nonlinear processes and 
show good agreements between, our calculated results 
and experimental observationscd as well as theoretical 
resultsED reported earlier. These experimental results are 
otherwise explainable only when many-electron calcula- 
tions are doneEEl. This is an indication that the effective 
medium theory is capable of representing the most im- 
portant many-electron effects in such systems in a direct 
and physically attractive way. Recent calculations by Das 
Sarma and HwangLj also show the importance of the bar- 
rier width in understanding the behavior of plasmons in 
MQWs, which is in agreement with our results. 

We take a MQW system as a layered medium with al- 
ternate layers of width Lw = ad and Lb = (1 — a)d 
where d is the period of the MQW and a is a frac- 
tion. For d <C A, the wavelength of light, the composite 
medium canJae described by an effective dielectric func- 
tion, iij(uj)cB. The calculation of the nonlinear response 
has to be performed in a self-consistent manner since the 
electric field in the well region, Ew and in the barrier 
region, Eb depend on the dielectric functions as 



Ew = e B E/ {e B a + e w (l - a)} . 



(12a) 



X (-w) = - 
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(m* aVw (E)f a (E) - m* bVw (E)f b (E)) 
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l+(hu- E b 



) 2 r 2 - 2 



■ Az ab Zba\Ew\ 2 (riT 2 



-dE (11) 



where a and b are two consecutive energy levels (b > a), 
Ti = h/T\, and r 2 = T1/T2 with T\, T 2 being the lifetime 
and the dephasing time, respectively. 



Electron-electron interaction 

It is generally accepted that the Coulomb interaction 
between the carriers plays an important role in the en- 
hancement of optical nonlinearity and the blue shift of 
intersubband resonances. For many systems, e.g. a 
metallic sphere, or a metal-insulator interface simple elec- 
tromagnetic theories of surface plasmons and surface 
polaritons rbave been important in understanding this 
interactionci The effective medium theory is a similar 
one which we will apply for MQWs in this section. The 
direct physical nature of this theory allows us to predict 
the plasmon frequency-dependence on physical param- 
eters, such as the width and the dielectric constant of 
the barrier material. For low carrier concentration the 
modification from single electron approximation is not 
so important, on the other hand, for high, carrie r con- 
centrations collective effects are importanttSllallj. We 
have already shown that the effective medium theooz. 
can give a good representation of thc-collective effectst£l 
in linear absorption spectra of ISBT&3. Here we extend 



Eb = e w E/ {e B a + ew{l ~ a)} , 



(12b) 



where E corresponds to the average electric field in the 
^-direction given by, 



E = aE w + (1 - a)E B . 



(13) 



The function ew, the dielectric permittivity of the well 
layer, in turn, depends on Ew according to ew = tbuik + 
^intersubband where e int er subband is calculated from Eq. 
(pd|). For simplicity we have assumed that eg, the di- 
electric permittivity of the barrier layer is a constant in 
the frequency range of interest. In case of z-polarized 
light, the electric displacement vector across the inter- 
face is continuous, i.e. 



ewEw — £bEb = £E. 



(14) 



From Eqs. fll2j), (|T3|), and ( O ) we obtain the average 
dielectric constant of the composite medium as 

tw^B 



From Eqs. 
relation: 



P 2 C 2 x 3 



aeB + (1 — a)ew 
l]), ([l4]), and ( [i"5| ) we obtain the following 



2Re (PCPoe^) + 2P 2 CB - yC 2 
B 2 P 2 + f3 2 + 2Re (BPfioe^) 
2BCy] x - yB 2 = 0, (16) 



where 
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P — £B + Efmifc + 1; 

1 — a 

B = l + (hLu- E ba fTz 2 : 

c = iz ab z ba (rir 2 ) _1 , 

e 2 

A)e^ = — (n sa ~ n sb ) \z ba \ 2 (Tilo - E ba - iT 2 ) T^ 2 , 



x = \E W \ 2 . 

In the above n s , is the sheet carrier concentration of the 
i-th energy level. Equation (|l6|) admits multivalued so- 
lutions for intensity in the well region as a function of 
the average intensity, \E\ 2 for several values of physical 
parameters. 

The value of the absorption co-efficient can be calcu- 
lated from the imaginary part of e where the solution 
of Eq. ([l6]) is used to calculate ew- The value of the 
nonlinear refractive index for large detuning or below the 
satuaration intensity can be calculated from 

X (-u;) = x« (- W ) + 3 X ?L(-u>)\E w \ a , (17) 

where xS (— u>) is the linear susceptibility due to ISBT. 
In this case, the collective effect can be intrpduced fol- 
lowing the formulation as reported elsewhereta. 

III. NONRESONANT xKL(O) 

{ 31 

The nonresonant third order nonlinearity Xzzzz (0) due 
to ISBT has been calculated using Eq. (||) for .aii 
InAs/AlSb MQW system used by War burton et alM 
The well and the barrier widths are 18 and 10 nm, respec- 
tively, and the carrier concentration is 2.61 x 10 12 cm~ 2 
at 10K whereby the second subband is populated. The 
values of the physical parameters for the calculations, 
and the results are given in Tables I and II, respectively. 

TABLE I. Physical parameters used in the calculations 

Material Systems mj^/mo m* B /mo aw olb Ecb — Ec 

(eV- 1 ) (eV" 1 ) (eV) 

InAs/AlSb 0.025 0.260 2.050 0.223 1.350 

GaAs/GaAlAs 0.067 0.0919 0.573 0.295 0.261 

toq is the free electron mass. 



TABLE II. Calculated values of x*iU(0) in e.s.u. 



Using Eq. (H 




Using Eq. (Jicj) 


Parabolic 


Nonparabolic 




-9.80 x lO"' 


-5.08 x 10"' 


-5.27 x 10" a 



The values of Xzzzz(0) are negative in all cases whereas 
the low-frequency mobile carrier contribution to x^ f° r 
bulk InAs is positive. It is found that Xzzzz(0) is non- 
zero even for the parabolic band approximation, which 
suggests that nonlinearity is mainly due to confinement 
of electron wave function within the quantum well region. 
When nonparabolicity is included, the dipole matrix el- 
ements do not change muchEH, but the zone centre ISBT 
energies decrease compared to those for the parabolic 
band approximation, and continue to decrese with in- 
crease in the in-plane momentum. As a result, both the 
terms within the parenthesis in Eq. (|^) increase. But, 
the rate of increase of the first term which is smaller of 
the two terms, is more than that of the second. Hence 
|xEL(0)| decreases. 

It is instructive to compare these results with the val- 
ues of Xzzzz (0) calculated for infinitely deep QWs using 

Eq. We see that the value of |xizL:(0)| is much 

larger than that obtained by using Eq. ( JlOj ) for infi- 
nite barrier model. This is expected for two reasons. 
Firstly, since a real-life QW system has a finite barrier 
height, the electron wavefunction penetrates into the bar- 
rier region. As the electron spends part of the time in 
the barrier layer, the effective well width increases, and 
|Xzzzz(0)| being proportional to the ninth power of the 
well width, for finite barrier heights, its value increases 
by an order of magnitude even for the parabolic band ap- 
proximation. Secondly, in the calculation using Eq. ( p^| ) 
partially occupied subbands are not included, i.e. only 
the transitions at or near the zone centre are considered. 
This precludes transitions from the first to the second 
subband when the iatter is populated. However, such 
transitions do occuiIIj for kf2 < k± < kfi where kfi and 
kf2 are the Fermi wave vectors for the first two subbands. 

Such transitions have been included in Eq. (O), which 

(3) 

leads to an increase in the value of |xzzzz(0)|- 

From these results we may conclude that nonlinearity 
due to ISBT is influenced by both the band nonparabol- 
icity as in bulk semiconductors, and the quantum con- 
finement of the electron wavefunction. The main contri- 
bution is from the quantum confinement which is sub- 
stantially modified by nonparabolicity. 

IV. THIRD HARMONIC GENERATION, xSL(3w) 

(3) 

The real and the imaginary parts of Xzzzz due to THG 
are shown in Fig 1. The calculation have been done for 
InAs/AlSb MQW system by sum-over -states method. 
There are 14 bound levels, all the levels are taken into 
account in calculations. It is seen that the value x^ (3w) 
converges after taking 10 bound levels. Because sum of 
the oscilator strengths for the bound states becomes close 
to unity, the effect of continuum may be safely neglected. 

(3) 

In Fig 2. | X zzzz(3w)| is shown, which is a measurable 
quantity. There are 8 peaks due to w, 2ui, and 3w res- 
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onances from the two occupied subbands. The u> , 2u>, 
and 3w resonances due to transitions from the ground 
subband occur at frui — 77, 95, 25, and 105 meV, respec- 
tively. The resonances due to transitions from the first 
excited level occur at hco = 110, 117, 36, and 120 meV, 
respectively. The peaks at 25 and 36 meV are very small, 
as shown in the figure. It is seen that around 117 meV, 
i.e., corresponding to the 10.6 /im region, a very high 
value, about 4.4 x 10~ 4 e.s.u. of |xSzz(3w)| is obtained, 
which is the largest value of nonlinearity due to THG ever 
reported for a single symmetric QW system. Since this 
MQW system already exists, it is worthwhile to do ex- 
periments to establish high nonlinearities of purely elec- 
tronic origin. We note that when depolarization effects 
are included, one-photon resonances would shift from the 
calculated positions as in one photon spectra. 



V. NONLINEAR REFRACTION AND 

(3) 




0.02 0.04 0.06 0.08 0.10 0.12 0.14 

Photon Energy (eV) 

FIG. 1. Dependence of x (3w) on photon energy for 
the InAs/AlSb MQW system. The dotted and the solid lines 
show the real and the imaginary parts, respectively. 




Photon Energy (eV) 

FIG. 2. Dependence of Ix'' 3 ' (3w) | on the photon energy 
for the InAs/AlSb MQW system. 



ABSORPTION, X zLz{ 



The real and the imaginary parts of Xzzzz{— w) are 
shown in Fig. 3. The imaginary part describes TPA. 
Here we obtain four major peaks due to transition from 
the ground and the first excited state. 




Wavelength (|im) 

FIG. 3. Dependence of \^\~Ld) on the incident wave- 
length for the InAs/AlSb MQW system. The dotted and the 
solid lines show the real and the imaginary parts, respectively. 

From the figure it is seen that at about 10.6 /jm there 
is a resonance due to two photon absorption between the 
second and the fourth subbands. The magnitude of third 
order optical nonlinearity corresponding to this TPA is 
also very high, of about 4 x 10 e.s.u. It is compara- 
ble to the largest reported nonlinearity at 10.6 /xm for 
Hgo.84Cdo.i6Te (Ref. 32). It should be noted that this 



resonant nonlinearity is of practical interest because no 
linear absorption occurs at this wavelength. For all op- 



tical switching devices the ratio, 



Mx (3) ) 

/m(x< 3 ') 



representing 



the phase of the nonlinearity is important, and should 
have a value greater than 2 (Ref. ^5 18) for device appli- 
cations. This ratio is plotted in Fig. 4, and we see that 
the condition is satisfied near the TPA wavelength. We 
have earlier predictedliS intrinsic optical bistability in this 
sample considering electron-electron interaction through 
effective medium theory. 

The real and the imaginary parts of Xzzzz(-oj) of 
GaAs/GaAlAs MQW system are plotted in Fig. 5. The 
physical parameters used in the calculations correspond 
to the sample used by Craig et alEd are listed in Table I. 

The results for the higher carrier concentration (N s — 
1.1 x 10 11 cm -2 ) are given in this figure. Since only one 
energy level is populated here, two peaks, one each for 
one photon and two photon resonances are observed. The 
high value of x (3) (-^) ( 2 - 2 xl(T 2 e.s.u.) at the TPA 
frequency is due to the large well width of 40 nm and 
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high carrier concentration. Since 



Mx (3) ) 

Jm(x< 3 >) 



> 2 (see Fig. 



6) near the TPA frequency, this particular QW system is 
also expected to be useful for optical switching. 




Wavelength (\im) 



Re ( X (3) (-^))//m (x (3) (^)) 
on the incident wavelength for the InAs/AlSb MQW system 



FIG. 4. Dependence of 




0.006 0.008 0.010 0.012 0.014 

Photon Energy (eV) 

FIG. 5. Dependence of X ( — a; ) on the photon energy for 
the GaAs/GaAlAs MQW system. The dotted and the solid 
lines show the real and the imaginary parts, respectively. 

Using Eq. (0) we have also calculated \xz)zz{— u>)\ fpii 
GaAs/GaAlAs MQW systems used by Walrod et al.El 
and obtained a set of values of 3.65 x 10~ 5 and 3.91 
x 10~ 5 e.s.u. which are in good agreement with their 
experimental results, 3.3 x 10 -5 and 4.25 x 10 -5 e.s.u., 
respectively. 




0.008 0.010 0.012 

Photon Energy (eV) 



FIG. 6. Dependence of Re ( X (3) (-uj)) /Im (x (3) (-w)) 
on the photon energy for the GaAs/GaAlAs MQW system. 



VI. INTENSITY-DEPENDENT ABSORPTION 

To calculate the intensity-dependent absorption 
consider the GaAs/AlGaAs system used by Craig et al.Eil 
The values of transition electric dipole moments and the 
dephasing relaxation time, are taken from the exper- 
imental workEl. We also assume T\ = T%. In the effec- 
tive medium theory the intensity in, and the intensity- 
dependent dielectric function of the well region depend 
on each other through Eqs. ( [Tl| ) and (|l2]). These were 
solved self-consistently at each intensity to obtain e whose 
imaginary part determines the absorption by the compos- 
ite medium. These results are shown in Fig. 7. 




Frequency (cm ') 

FIG. 7. Dependence of scaled absorption on the frequency 
of the incident radiation for the GaAs/GaAlAs MQW system 
for Ti = Ta. The intensities of the radiation for the highest 
to the lowest peaks are 10" 3 , 100, 200, 500, 800, and 1000 
Wcm -2 , respectively. 



Here, the absorption co-efficient, multiplied by a scale 
factor, hnoc/q where no c, and q are the refractive index 
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of bulk GaAs, velocity of light, and the electronic charge, 
respectively is plotted in the ordinate. The asymmetry in 
the lineshapes of the absorption curves is known to occur 
due to the collective effectO and is well reproduced in our 
calculation. We note that the observed position of the 
linear absorption peak at 81 cm -1 is blue shifted from 
the calculated one electron peak position at 72 enr 1 , 
which can be explained as the depolarization shift with 
the help of the many electron theory. This shift has also 
been accurately reproduced in our calculations. Using an 
intensity-dependent T\ it is possible to get much better 
agreement with the experimental intensity-dependent ab- 
sorption spectrum for the sample having N s = 1.1 x 10 11 
cm~ 2 . For the other samples, a change in the absorp- 
tion frequency by 5 cm -1 is observed for high incident 
intensities. The discrepancy is due to the asymmetry 
in the QW shape caused by an applied DC electric field 
which has not been included in our calculations. We have 
also obtained similar results for the GaAs/GaAlAs sys- 
tem discussed by ZaluznyO. 

VII. CONCLUSIONS 

In this paper we have presented detailed calculations 
of the third order optical nonlinearity due to ISBT in 
semiconductor MQW systems. We first calculate the 
single particle linear and nonlinear response functions. 
The most important effects of electron-electron interac- 
tions i.e. the screening due to depolarization terms are 
included then through an effective medium theory. In 
calculating the single particle response all the essential 
complexities, such as the finite barrier height and the en- 
ergy band nonparabolicity are included. Very high values 
of x««(-w) (4 x 1CT 4 e.s.u.) and X ^L(3lu) (4.4 x 1CT 4 
e.s.u.) are predicted near 10.6 /xm for the InAs/AlSb 
MQW system. We find that in the low-frequency regime 
the nonlinearity is mainly due to quantum confinement 
of electrons - the effect of nonparabolicity being a modi- 
fication of this nonlinearity. 

The inclusion of electron-electron screening terms 
through an effective medium theory gives rise to several 
new features. In linear spectra the dominant absorption 
peaks are attributed to intersubband plasmons. The blue 
and the red shifts observed in the linear and the non- 
linear optical absorptions, respectively due to the many 
electron effect can be explained in terms of the effective 
medium theory. The asymmetry of the lineshapes of the 
nonlinear absorption spectra for intensities higher than 
the saturation intensity is also obtained. 
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